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Abstract 

We study the charged and neutral current weak interaction rates relevant for 
the determination of neutrino opacities in dense matter found in supernovae 
and neutron stars. We establish an efficient formalism for calculating differ- 
ential cross sections and mean free paths for interacting, asymmetric nuclear 
matter at arbitrary degeneracy. The formalism is valid for both charged and 
neutral current reactions. Strong interaction corrections are incorporated 
through the in-medium single particle energies at the relevant density and 
temperature. The effects of strong interactions on the weak interaction rates 
are investigated using both potential and effective field-theoretical models of 
matter. We investigate the relative importance of charged and neutral cur- 
rents for different astrophysical situations, and also examine the influence of 
strangeness-bearing hyperons. Our findings show that the mean free paths 
are significantly altered by the effects of strong interactions and the multi- 
component nature of dense matter. The opacities are then discussed in the 
context of the evolution of the core of a protoneutron star. 



PACS number(s): 13.15.+g, 26.60.+C, 97.60Jd 
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I. INTRODUCTION 



The transport of neutrinos is an essential aspect of simulations of gravitational collapse, 
supernovae, protoneutron stars and binary mergers of compact objects. The neutrinos of 
all flavors emitted from newly formed neutron stars in supernova explosions |J are the 
only direct probe of the mechanism of supernovae and the structure of protoneutron stars. 
The most important ingredient of neutrino transport calculations in these simulations is 
the neutrino opacity at supra-nuclear density [0-0]. Yet, to date, calculations of neutrino 
opacities in dense matter have received relatively little attention compared to other physical 
inputs such as the equation of state (EOS). 

Both charged current absorption and neutral current scattering reactions are important 
sources of opacity. Neutral current processes involve all flavors of neutrinos scattering on 
baryons and leptons. While scattering from electrons is important for energy and momen- 
tum transfer in the process of thermalizing neutrinos, for both energy and lepton number 
transport, neutrino-baryon scattering and absorption are the dominant processes. 

Earlier work on neutrino interactions in matter at supra-nuclear densities that have 
shaped our discussion in this paper are in Refs. [9-32]. Various approximations have al- 
most always been made in the calculations of neutrino cross sections. These approxima- 
tions concern the degrees of degeneracy or relativity, the composition, or the effects of the 
baryon-baryon interactions. The cross section for neutrino-nucleon interactions has only 
been calculated exactly for noninteracting gases by Schinder [^J and for neutrino-electron 
scattering by Mezzacappa and Bruenn || . Other workers have developed limiting expres- 
sions for noninteracting gases for the purely elastic case, the completely degenerate case, 
or the completely nondegenerate case. In fact, most numerical simulations of supernovae, 
protoneutron star evolution, and binary neutron star coalescence, have employed limiting 
expressions, derived from those for noninteracting nucleonic matter, in which interactions 
were included by simple scaling factors |2|,(|,[7]|| . The results may be inconsistent with the 
underlying nuclear matter EOS. 

Only a few attempts exist [ J^C , 17 , IS , 27 , 2£ ] in which the effects of strong interactions on 
neutrino opacities were considered. But these studies were not performed for the interme- 
diate degeneracies that are often encountered in astrophysical environments. In the case 
of neutral current opacities, the effects of stong interactions were investigated for nonde- 
generate nuclear matter by Sawyer JT(| and for degenerate neutron matter by Iwamoto & 



Pethick |17|]. Both of these efforts treated nucleons in the nonrelativistic limit and predicted 
increases in the mean free paths by factors of ~ 2 — 3, for densities in the range of 2 — 4 
times the nuclear saturation density (n = 0.16 fm~ 3 ). Subsequently, relativistic calculations 
based on effective Lagrangian models for hot, but neutrino-poor neutron star matter were 
performed by Horowitz and Wehrberger ||27|| . The neutral current differential cross sections 
were calculated using linear response theory and mean free paths were found to be factors of 
1.5 — 2 times that for noninteracting nucleons. In the case of charged current interactions, 
the effects of interactions have been considered by Sawyer [|HJ for nondegenerate matter and 
by Goodwin and Pethick [0 for degenerate matter. 

Schinder's exact results for noninteracting nucleons are expressed in terms of a lengthy 
series of Fermi integrals. We have found, however, a simplification to Schinder's results that 
not only expresses them in numerically simpler fashion as Polylogarithmic functions, but are 
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also easily generalized to the case of interacting matter. In contrast to previous work, we 
calculate both neutral and charged current opacities including effects of interactions arising 
from the underlying EOS. 

Another aspect of opacities which recent work has emphasized concerns the importance 
of including the multi-component nature of dense matter on neutrino opacities. Prakash 



et. al. |28] , and Reddy and Prakash PH,|33| have identified neutrino-hyperon absorption and 
scattering reactions as being important new sources of opacity. These include absorption 
involving the A and S~ hyperons and scattering involving the £~ hyperon. These could 
play important roles in calculations of the neutrino signature of an evolving protoneutron 
star with hyperons ||. 

In this paper, we perform neutrino opacity calculations for interacting matter of ar- 
bitrary degeneracy and composition at supra-nuclear densities. Interactions between the 
leptons are negligible and will be ignored. Strong interactions between the baryons, how- 
ever, significantly alter the energy spectra from their noninteracting forms, especially at 
high density. One of our objectives is to explore the extent to which interactions among the 
baryons affect the neutrino cross sections. We will separately consider potential models that 
are inherently nonrelativistic and relativistic field-theoretical models. Appendices A and B 
consider the formulation of these models in some detail. We also include effects due to the 
multicomponent nature of interacting matter and the possible presence of strangeness in 
the form of hyperons. Where possible, we provide analytical expressions for both charged 
and neutral current differential scattering cross sections for a given EOS at a fixed density, 
temperature and lepton concentration for given incoming neutrino energy and scattering 
angle. The most detailed transport codes 0, which solve the full Boltzmann transport 
equation, require differential cross sections. However, simpler transport methods ||-|| only 
need angle and/or energy averaged opacities, which can be usually expressed in the form 
of mean free paths. In this paper, to facilitate comparisons between our results and those 
involving various approximations, we present numerical results for mean free paths, leav- 
ing for subsequent publications detailed results for the general opacities. In particular, we 
present results for mean free paths for thermodynamic conditions relevant to the evolution 
of protoneutron stars. 

The paper is organized as follows. In Sec. II, we summarize the basic relations needed to 
calculate both absorption and scattering opacities in hot and dense matter. We will consider 
in Sec. Ill the idealized situation in which the baryons are treated as nonrelativistic, the 
leptons are assumed to be massless, and the baryon-baryon interactions are neglected. In 
this case, the baryon energy spectra has the form E = p 2 /2M. Our results lead naturally 
to previously obtained limiting cases such as the degenerate, nondegenerate, and elastic 
approximations. We will next consider, in Sec. IV, the nonrelativistic Skyrme-like potential 
model HU, in which the spectrum is given by E = p 2 /2M* + Unr, where Unr denotes 



the density-dependent potential contribution and M* is the Landau effective mass, which 
is also generally density dependent. Next, in Sec. V, we will consider the effective field- 
theoretical Walecka-type model PSJ at the mean-field level in which the spectrum is E — 



p 2 + M* +Ur, where M*, the Dirac effective mass, and Ur, which accounts for interactions 
of the fields, are both generally density dependent. The kinetic parts in both nonrelativistic 
and relativistic approaches contain the effects of interactions insofar as they depend on 
density-dependent effective masses; further, in both cases, the momentum dependence of the 
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kinetic energy is formally identical to their noninteracting counterparts. The evaluation of 
the cross section is therefore similar to the case of noninteracting baryons for these particular 
models. The numerical results, however, are sensitive to the presence of interactions. In 
Sec. VI, we show how additional baryonic components, such as hyperons, affect the neutrino 
opacities. In Sec. VII, we compare our results with those of previous workers. In Sec. VIII, 
we apply our results to the neutrino opacities encountered in a particular astrophysical event: 
that of a deleptonizing and cooling protoneutron star. Finally, we summarize and provide 
an outlook in Sec. IX. Appendices A and B contain supporting material about the potential 
and field-theoretical models of dense matter. 



II. NEUTRINO CROSS SECTIONS 

The total opacity of dense matter to neutrinos has contributions from both neutral and 
charged weak currents. Neutral currents are involved in neutrino-baryon scattering while 
charged currents are involved in neutrino-baryon absorption reactions. We have considered 
neutral currents in Ref. [33|, in which convenient expressions for the differential and total 



scattering cross sections were established. Here, we will concentrate upon deriving similar 
relations for the charged current absorption reactions. Further, the formalism developed 
here may be easily generalized to include neutral current reactions so that the present work 



in effect extends and replaces Ref. f33j. 



The neutrino energies of interest to us are less than a few hundred MeV; we may thus 
write the relevant interaction Lagrangian from Wienberg-Salam theory [|3B} -|3~8|1 in terms of 
a current-current interaction: 

£«=^ W for u l + B 2 ^l + B 4 (1) 

<^nt=% W fOT H + + (2) 

where Gf — 1.436 x 10~ 49 erg cm -3 is the Fermi weak coupling constant and the Cabibbo 
factor C = cos 6 C for change of strangeness AS" = and C = sin 6 C for AS = 1 . The lepton 
and baryon weak charged currents are: 

= (! - 75) , 3w = i J 4^(9v-9Al5)'^2- (3) 
Similarly, the baryon neutral currents are given by 

^ = VV7 M (! - 75) A , 3z = i^il 1 * i c v - c a7s) ^2 , (4) 

where 2 and 4 are the baryon initial state and final state labels, respectively (these are 
identical for neutral current reactions). Other particle labels and four- momenta are as 
shown in Fig. 1(a) for the charged current reaction and Fig. 1(b) for the neutral current 
reaction. The vector and axial-vector coupling constants gy and g& are listed in Table I 
for the various charged current reactions of interest. Similarly, the couplings cy and ca for 
the neutral current reactions are listed in Table II. Generally, the \x and r neutrino charged 
current reactions are kinematically suppressed; \x and r neutrinos are thermally produced so 
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that their energies are of order T « m M . On the other hand, neutral current reactions are 
common to all neutrino species and the neutrino-baryon coupling is independent of neutrino 
flavor. Neutrino coupling to leptons in the same family is modified since the scattering may 
proceed due to both W and Z exchange; the couplings shown in Tables I and II reflect this 
fact. Numerical values of the parameters that best fit data on charged current semi-leptonic 
decays of hyperons are f39|: D=0.756, F=0.477, sin 2 W =O.23 and sinfl c = 0.231. Note that 
the amplitude for the strangeness changing charged current is suppressed by the factor sin 9 C . 
These couplings follow from SU(3) flavor symmetry for the baryons and the quark model. 
Corrections arising due to explicit SU(3) breaking terms have been recently investigated HO 
and in some cases are about 10-30%. 

The cross section per unit volume of matter (or equivalently the inverse collision mean 
free path) may be derived from Fermi's golden rule and is given by 



V 



d 3 p 2 f d 3 p 3 f d 3 p. 



(2tt) 3 J (2tt) 3 7 (2tt) 3 

x / 2 (£ 2 )(1 



(2tt) 4 5 4 (P 1 + P 2 - P 3 - P 4 ) W fi 
f 3 (E 3 ))(l - f 4 (E 4 )) , 



(5) 



where Pi = (Ei,pi) denotes the four-momentum of particle i (particle labels are as shown in 
Fig. 1(a)) and the transition rate Wfi is 



W 



ft 



< \M\ 2 > 
2 4 EiE 2 E 3 E 4 



(6) 



Above, \A4\ 2 is the squared matrix element and the symbol < • > denotes a sum over final 
spins and an average over the initial spins. A common expression for both scattering and 
absorption may be written: 



W 



fi 



(V + A) 2 {1 -v 2 cos 6>i 2 ) (1 - v 4 cos9 34 ) 



+ (V--4) 2 (1-V 2 COS#23)(1-W4COS0 



14 



- (V 2 -^4 2 



| M 2 { 
E2E4 



cos 6 



13; 



(7) 



where the vector and axial couplings V and A, in the case of absorption, stand for Ggy and 
CgA, respectively. For the reactions of interest, g v and g& are listed in Table I. Similarly, for 
the scattering reactions of interest, V and A stand for cy/2 and c^/2, respectively, which are 
listed in Table II. The particle velocities are denoted by Vi = Pi/Ei, and the angle between 
the momentum vectors pi and p~j is denoted by 9ij. Further, M is the bare nucleon mass. 
The functions fi(E{) in Eq. @ denote the particle distribution functions, which in thermal 
equilibrium are given by the Fermi-Dirac functions 



fi(Ei 



1 + exp 



E; 



T 



(8) 



where E^ are the single particle energies, /ij are the corresponding chemical potentials, and 
T is the temperature. 

In general, the single particle energies and chemical potentials depend on the ambient 
matter conditions, i.e., the density and temperature, and also on the interactions among the 
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various particles. The various chemical potentials are determined by the conditions of charge 
neutrality and, in all but the most extremely dynamical situations, chemical equilibrium. 
In some astrophysical situations, such as in the late stages of core collapse and during the 
early stages of the evolution of a protoneutron star, neutrinos are trapped on dynamical 



times within the matter [f£ip2| and chemical equilibrium is established among the baryons 



and leptons. In this case, the chemical potentials satisfy the relation 

VB 2 ~ VB A = He ~ Hv B - (9) 

These situations are characterized by a trapped lepton fraction Yl = Y e + Y Ue , where 
Y e = (n e — n e +)/nB and Y Ve = (n Ue — n Ve )/nB are the net electron and neutrino fractions, 
respectively. The evolution of a protoneutron star begins from a neutrino-trapped situation 
with Yl w 0.4 to one in which the net neutrino fraction vanishes and chemical equilibrium 
without neutrinos is established. In this case, the chemical equilibrium is modified by setting 
fj, Ve = 0. In all cases, the condition of charge neutrality requires that 

i i 

where the superscript's (±) on the number densities n signify positive or negative charge. 

Although neutrino opacities are required for a wide range of densities, temperatures, and 
compositions, for the most part we will display results for two limiting situations, namely 
beta equilibrium matter with either Yl = 0.4 or Y v = 0. These are situations encountered 
in the evolution of a protoneutron star K3|, as discussed further in Sec. VII. 



III. NONRELATIVISTIC NONINTERACTING BARYONS 

For baryon densities % < 5no, where no = 0.16 fm -3 is the empirical nuclear equilibrium 
density, and in the absence of interactions which could significantly alter their effective 
masses, baryons may be considered as nonrelativistic. The expression for W/i in Eq. (^) 
then simplifies considerably, since the baryon velocities < 1. In this case, the terms 
involving the baryon velocities may be safely neglected. However, the term involving the 
angle between the initial and final leptons remains. For reactions involving nucleons, this 
term gives a small contribution, since it is proportional to V 2 — A 2 . For simplicity, and to 
make an apposite comparison with earlier results in which this term was also neglected, we 
drop this term in this section, but will return to a more complete analysis in the succeeding 
sections. 

Under these conditions, the transition rate Wfi becomes a constant, 

W fi = G 2 F (V 2 + 3A 2 ), (11) 

independent of the momenta of the participating particles, and the differential cross section 
is given by 
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where the three-momentum transfer q = p[ — p~* 3 , so that q = \q\, and the energy transfer 
go = E\ — E 3 . The function S(qo, q), the so-called dynamic form factor or structure function, 
characterizes the isospin response of the (nonrelativistic) system. It is simply the total phase 
space available to transfer energy q and momentum q to the baryons. We note that the 
differential cross section is needed in multi-energy group neutrino transport codes. However, 
more approximate neutrino transport algorithms often only require the total cross section 
as a function of the neutrino energy. The cross section per unit volume given in Eq. (^J) then 
simplifies to 

^ = G 2 F (V 2 + 3A>) J |^ (1 - fsm)S(q , q), (13) 

s{q(h q) = 2 I^I§^ (27r) V(Pl + P2 ~ Ps - w^X 1 - M E *» • ( 14 ) 

The total cross section given by Eq. ( |I3"D can be recast as a double integral in (q , q) space 
using d 3 p 3 = 2irq(E 3 /Ei) dq dq. Since E 3 ranges between and 00, the limits of q are —00 
and Ei. The limits of q are obtained by inspecting the relation q 2 = E 2 + E 2 — 2EiE 3 cos 9 
for cos #13 = ±1. Thus, |go| < q < 2Ei — q . One finds 
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v !t2 (V 2 + 3^1 2 ) J_y q o^(l-f 3 (Es))J M °dqqS(q ,q). (15) 

This, of course, applies to both scattering and absorption with appropriate changes of par- 
ticle labels and coupling constants. 

The integrals in Eq. fll4]) can be performed analytically and the result expressed in closed 
form for the noninteracting case and for certain models of interacting matter. The integral 
over the final state momentum in Eq. (Jl4j) may be performed by exploiting the momentum 
delta function to obtain 

S(q ,q) = ^ fd 3 p 2 5(q + E 2 -E,)f 2 (E 2 )(l-UE 4 )). (16) 

We note that E 4 = (p 2 + q) 2 /2M, and we ignore the difference between M 2 and M4 as it is 
small compared to other energy scales at high density for noninteracting matter. We may 
rewrite the energy delta function in terms of the angle between p 2 and q: 

5(q + E 2 -E i ) = — 5(cos 9 - cos 9 )e(p 2 2 - p 2 _) , (17) 
p 2 q 

where 

a M ( q 2 \ 2 M 2 ( q 2 \ 2 

COSe0 = p^q{ q0 ~2M)' P - = ^{ q °-2M) (18) 

and 6(x) = 1 for x > and zero otherwise. Substituting these results in Eq. fll6|) and 
performing the angular integrals we obtain 



S(q , q) = — r dp 2 p 2 f 2 (E 2 ) (1 - f A (E 4 )) . (19) 
Tiq Jp- 
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The remaining p 2 integral is performed by using the relation 



dx 



1 



In 



1 + exp(.x) 



J l + exp(x) 1 + exp(— x — z) 1 — exp(— z) 1 + exp(x + z) 

Thus, the free gas isospin density response function is given by 

M 2 T 



(20) 



S(qo,q) 



irq 



1 — exp(— z) 



(21) 



where 



Z = T » ^ — ^2 — ^4, 



= In 



1 + exp((e_ - fi 2 )/T) 
1 + exp((e_ + g - ^)/T) 



e_ 



p- i (go - g 2 /2M) ; 



2M 



q 2 /2M 



(22) 



This result generalizes a result obtained earlier for noninteracting symmetric nuclear mat- 
ter in which /t = 0, to the case of asymmetric nuclear matter for conditions of arbitrary 
degeneracy. This result, which we further generalize to include nuclear interactions in the 
next section, proves to be the key to being able to efficiently calculate opacities. These 
results are easily specialized to the case of scattering by noting that particle labels 2 and 4 
are identical; thus fa = fa- The integrals in Eq. ( fCEf) , even with an analytical expression 
for S(qo,q), require numerical evaluation; closed form expressions for arbitrary degeneracy 
cannot be obtained. However, in some limiting cases these integrals become analytic and 
correspond to results obtained earlier and which are often used in astrophysical simulations. 

Highly Degenerate Baryons (fa/T ^> 1): In this situation, the participating particles all lie 
close to their respective Fermi surfaces. In this case, the q integration may be performed 
trivially, since the factor (1 + may be replaced by Q(fa — e_) or, equivalently, Q(q — 

(Pf 2 ~ PfJ), where B(x) = 1 for x > and zero otherwise. The integral to be performed is 
then 



2E 1 -q 



1o\ 



dq Q(q - (p F2 -p F J) 



IAI 



dq Q(q- (pf 2 ~ PfJ) 



(23) 



where in writing the rightmost relation, we have set go — — /*> since the exponentials in the 
g integral are highly peaked at this value at low temperature. Thus, 



2E 1 Q(fi - (p F2 - pf a )) : for p, > p F2 
(2Ei - fx - pf 2 + Pf 4 )©(2-Ei - fi - Pf 2 + Pf a ) ■ for fi < Pf 2 



Pf 4 
Pf 4 



(24) 



The upper limit on the remaining g integral can be replaced by +oo since the integrand 
vanishes exponentially for positive values of go due to final state Pauli blocking of the electron 
degeneracy. With this substitution, it is straighforward to perform the go integral by noting 
that 



dz 



1 



1 — exp (—z) 1 + exp (z + rj) 



tc 2 + rf 



1 + exp (77) 



(25) 
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The final result for the cross section per unit volume in the degenerate approximation is 
given by 



G 



F (V 2 + 3A 2 )M 2 T 2 E (Ex + fi) 



cr(E 1 ) 

V 4vr 3 

with S = Q(pf 4 + Pf 3 ~ Pf 2 ~ Pf 1 
Pf 4 + Pf s - Pf 2 + Pf 1 



2 I Ei — fa 



1 



l + exp((fn- Ei)/T) 



+ 



2Ei 



©(PFj - \pf 4 +Pf 3 -Pf 2 \)- 



(26) 



In the above expression all terms proportional to T j Hi are neglected as they are small. Fur- 
ther, if we were to assume that the incoming neutrino energy Ei were equal to the neutrino 
chemical potential, the factor (E\ + ft) could be replaced by He- With this substitution, the 
above result coincides exactly with that derived earlier by Sawyer and Soni [|T5|]. 
In the case of scattering, this expression simplifies to 



°{Ei 
V 



G 



167T 3 



' 'c^ + 3c\)M 2 T 2 Ei 



i Ei - hi 



l + ex P ((Hi-Ei)/T) 



(27) 



An analogous result was derived earlier by Iwamoto and Pethick |Tj 

Note that the energy e_ arises due to the kinematical condition that ensures three- 
momentum conservation. For e_ ^> H2, in the degenerate limit, the phase space rapidly 
vanishes. Thus, at low temperatures, this leads to the condition 



q > \pf 2 -PfA 



(28) 



which is the threshold condition for the so-called direct Urea process [}45|J . Note that the 
maximum possible momentum transer is q — E\ +Pf 3 - For a free gas in beta equilibrium 
with zero neutrino chemical potential, the condition in Eq. (53) is usually not satisfied, since 
the neutron Fermi momenta are usually much larger than those of protons and electrons. 
In contrast, in the case when neutrinos are trapped, significantly larger proton fractions are 
favored in beta equilibrium condition. This enables the threshold condition to be easily ful- 
filled. Similarly, strong interactions also tend to increase the proton fraction, which in some 
cases allows Eq. (|28| ) to be satisfied even in the vicinity of nuclear densities The factor 
1 + z naturally accounts for the threshold- like behavior with decreasing temperature. 



Nondegenerate matter (Hi/T —1): In the nondegenerate limit, one has 

>2 - e_ 




cxp 



[1 - exp(-z)] 



(29) 



after expanding the logarithmic terms in Eq. (22) to leading order. Thus, S(qo,q) ~ 
(M 2 T/7rg) exp((/i2 — e_)/T). Assuming further that the effects due to final state Pauli 
blocking may be neglected and the relevant energy transfer is small, go *C Ei, one finds 



a{Ei) G 2 F (V 2 + 3A 2 )M 2 T //x 2 \ r^i 



V 



47T 3 



vx])[ T)Lj q ° C dqexp Ht 



(30) 



In this case, the go integration may be performed first. To leading order in q and T/M the 
result is g(27rT/M) 1 / 2 . The remaining q integration is elementary and we obtain 
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cr(£i) 
V 



G 2 F 

7T 



(V 2 + 2>A 2 )El n 2 , 



(31) 



where n 2 = 2(MT/2n) 3 / 2 exp(fx 2 /T) is the neutron number density in the nondegenerate 
limit. In this limiting situation, the total cross section is simply that on a single baryon 
times the baryon number density M. 



The Elastic Approximation: We can also derive the result of the so-called elastic approxi- 
mation || , in which it is assumed there is no energy or momentum transfer to the nucleons. 
The effects of final state Pauli blocking and nucleon degeneracy are still fully included. In 
this case, the cross section may be obtained by considering the response function for q — > 
and q — > 0: 



S(q ^0,q^0) 



2ir8(E 1 
2ttc%o)- 



En 



d 3 p 2 
(2tt) 3 
n 2 - n 4 



f 2 (E 2 )(l-f 4 (E 2 )) 



exp((/i 4 - fi 2 )/T) 



(32) 



where n 2 and n 4 are the number densities of species 2 and 4, respectively. In obtaining this 
result, the identity 



f 2 {E 2 ){l-U{E 2 )) 



f2{E 2 )-h{E 2 ) 
1 - exp((/i 4 - n 2 )/T) 



was used. One then finds 



n 2 - n 4 



V 



7T 



1 - exp((/i 4 - n 2 )/T) 



El (1 - h{E x )) . 



(33) 



(34) 



To specialize this result to the case of scattering, we need to take the limit (fi 2 — /z 4 ) — » 0: 



V 



^(4 + 3c 



dn-. 



El (1 - / 3 (£i)) • 



(35) 



We now present some results for noninteracting matter containing nucleons, electrons, 
and trapped neutrinos for ub and T chosen to lie in the range (0.25 — 5)uq and (5 — 60) MeV, 
respectively. For the beginning of the evolution of a nascent neutron star, it is appropriate to 
consider the neutrino trapped case in which the electron lepton fraction Yl = Y e + Y Ue is held 
fixed at 0.4. At later times the neutrino chemical potential vanishes and the matter contains 
thermally produced pairs; under these circumstances, Y Ue and fi Ue vanish and Y L = Y e is 
determined by charge and chemical equilibrium. 

Fig. 2 shows the composition and the chemical potentials versus the density (in units 
of the nuclear equilibrium density no) for T = 5, 30, and 60 MeV for the neutrino-trapped 
case. The major effect of trapping is to keep the electron concentration high so that matter 
is more proton rich in comparison to the case in which neutrinos do not carry net lepton 
number; thus the threshold conditon for the charged current rate is easily satisfied for all 
relevant densities and temperatures of interest. There is little variation in the individual 
concentrations with either density or temperature. The individual chemical potentials and 
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A = Hn — Hp increase significantly with density, but are relatively insensitive to variations in 
temperature as leptons are degenerate except at the lowest densities. 

In the neutrino-trapped case, in which the neutrinos are degenerate, Pauli blocking en- 
sures that the transport processes involve only neutrinos close to the Fermi surface. Thus 
it is appropriate to calculate the neutrino mean free path for a neutrino energy given by 
E\ = Hu e - The neutrino absorption mean free path, A = (a/V)" 1 , is shown in Fig. 3, as 
solid lines for temperatures T = 5, 30, and 60 MeV as a function of density. As expected, 
the mean free path decreases with increasing density since /i e increases with density. For 
comparison, we also show the results of the three limiting expressions for the extremely de- 
generate (Eq. (P^), dashed curve), nondegenerate (Eq. (|3T|), dot-dashed curve), and elastic 
(Eq. (|34f) , long dashed curve) approximations. It is evident that, for the densities and tem- 
peratures of interest to us, it is essential to include Pauli blocking of the final state electrons 
to obtain reliable results. This is chiefly due to the fact that electrons are degenerate. When 
all participating particles are extremely degenerate, Eq. (|2T)| ) provides an excellent approxi- 
mation to the exact result. The failure of the elastic approximation in this regime is related 
to the fact that the favored energy transfer is not zero as is assumed, but go — — Aj which is 
generally large in magnitude. At high temperatures and relatively low density, in which the 
nucleons are not extremely degenerate, the cross section in Eq. (0), based on the elastic 
approximation, reasonably reproduces the exact results. As expected, the nondegenerate 
approximation is only realistic for extremely low densities. 

Similarly, we have calculated the mean free path due to scattering reactions assuming 
Ei = fi Ue , and these are compared with the charged current results in Fig. 4 for T = 10 MeV 
as a function of density. The total scattering opacity, which has contributions from neutrons 
(dashed curve), protons (dot-dashed curve), and electrons (dot-long dashed curve), is shown 
as the long-dashed curve labelled A5. The results for electron scattering are actually taken 
from the relativistic formulation presented in Sec. V; they cannot be evaluated from the 
nonrelativistic formulation in this section because the electrons are highly relativistic. In the 
right panel, the ratio of scattering mean free path to that of the absorption mean free path 
is plotted. The solid curve is the ratio of the total scattering mean free path to that of the 
absorption mean free path, the dashed curve is the ratio of the mean free due to the reaction 
u e + n —>■ v e + n to that of absorption. u e + n — > e~ +p. Clearly, the scattering mean free path 
is dominated by neutrons since the electron and proton fractions are small. However, the 
net electron-neutrino mean free path in the trapped regime is dominated by the absorption 
reaction, for two reasons. First, the charged current matrix element is four times larger and 
second, the rate is proportional to /x e while the scattering rate is proportional to fj, Ve . 

In the neutrino-free case (Y Ue = 0), neutrinos still exist, but as pairs and do not carry net 
lepton number. Thus the reference to this phase as neutrino-free is somewhat misleading. 
Fig. 5 shows the composition and the chemical potentials. Compared to the neutrino-trapped 
case, Fig. 2, matter contains a smaller proton fraction. Consequently, the difference between 
the neutron Fermi momentum and proton (and electron) Fermi momenta are large, making 
simultaneous energy and momentum conservation impossibile. Thus, at low temperatures 
for all densities shown, neutrino absorption due the reaction u + n — > e~ + p is kinematically 
suppressed and is unimportant; the dominant absorption rates are those due to the modified 
Urea reaction v + n + X — > p + X + e~ where X is a bystander particle. The modified Urea 
reaction, however, gives a small contribution to the total opacity when compared with the 
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total scattering opacity and can be neglected. Although the direct reaction (u + n — > e~ +p) 
is kinematically forbidden at zero temperature, it exhibits a strong temperature dependance. 
We also note, as can be seen in Fig. 5, that beta-equilibrium and charge neutrality favors a 
larger proton fraction at higher temperatures, which enhances the importance of the direct 
process. In this sense, finite temperature effects on the neutrino absorption reaction are 
similar to those encountered due to the presence of a finite neutrino chemical potential. It 
increases both the neutrino energy and the proton fraction so as to simultaneously fulfill 
momentum and energy conservation. 

In Fig. 6, we plot the neutrino mean free path for neutrino-free matter with Y Ue = 0. 
Because the neutrinos are thermal, we choose E\ = 3T as the relevant energy for calcu- 
lating the mean free path. At low temperatures, the neutrino absorption mean free paths 
are very large and increase with density as matter becomes increasingly degenerate. For 
higher temperatures, this behaviour is not seen and the larger neutrino energy and a higher 
proton fraction both ensure that momentum and energy conservation for the direct reaction 
are easily satisfied. The extent to which the elastic and nondegenerate approximations re- 
produce the exact numerical results may be gauged from Fig. 6. Note that the degenerate 
approximation is not an appropriate limit because when the theta functions in Eq. ( f26|) are 
zero causing the mean free paths to be infinite; thus the degenerate limit is not plotted in 
this case. 

From Fig. 6 we see that the absorption rate shows a strong temperature dependance. In 
Fig. 7, we show the temperature dependance of the absorption and scattering mean free paths 
at nuclear saturation density. The individual scattering contributions are shown in the left 
panel together with the total scattering mean free path, which is shown as the long dashed 
curve. For this density, the scattering opacity is dominated by neutron scattering. The 
right panels show the ratio of the total scattering mean free path to that of the absorption 
reaction u + n — > e~ + p. For low temperatures, where the absorption is kinematically 
suppressed, scattering dominates the net opacity. However, with increasing temperature, 
absorption increases and eventually dominates for T > 10 MeV. 



IV. NONRELATIVISTIC INTERACTING BARYONS 

We turn now to the effects of strong interactions on the charged and neutral current 
neutrino cross sections. To begin, we will continue to use the nonrelativistic approxima- 
tion for baryons. Nonrelativistic potential model descriptions of dense matter are based 
on a two-body potential fitted to nucleon-nucleon scattering and a three-body term whose 
form is suggested by theory and whose parameters are determined by the binding of few- 
body nuclei and the saturation properties of nuclear matter (see for example, Ref. [§j| ). 
However, microscopic calculations of this type which encompass asymmetric matter at an 
entropy/baryon in the range s = 1 — 2 are not yet available. We therefore use a schematic 
potential model [47], which is designed to reproduce the results of more microscopic calcula- 



tions of both symmetric and neutron matter at zero temperature and which can be extended 
to asymmetric matter at finite temperature. In addition to employing local contact interac- 
tions, momentum dependent interactions arising from finite range exchange forces can also 
be considered in such an approach. The relevant details are given in Appendix A. 
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The calculation of the neutrino opacities is greatly simplified if the finite range interac- 
tions are approximated by effective local interactions. A notable feature of this approxima- 
tion is that the particle effective masses generally drop more rapidly with increasing density 
compared to the case in which the full nonlocality is retained. In this case, one retains only 
a quadratic momentum dependence in the single particle spectrum, so that it takes a form 
(cf. Appendix A) closely resembling that of a free gas: 



Ei{pi 



pj 

2M* 



n,p . 



(36) 



Both the single particle potentials Ui and the (Landau) effective masses M* are density 
dependent. Because the functional dependence of the spectra on the momenta is identical to 
that of the noninteracting case, the methods employed in the previous section are amenable 
to the incorporation of the effects of strong interactions. However, it is important now to 
retain the distinction between the masses of the nucleons, in particular the effective masses, 
M* and M*. 

The dynamic form factor in Eq. (|T4|), with the single particle spectrum in Eq. (|36"D, 
may be expressed in essentially the same functional form as in Eq. (|2l|) . The energy delta 
function can be recast in terms of the angle between p 2 and q: 



where 



5(q + E 2 -E 4 



, Ml 
COS#o = 

P2q 



Ml 

-^<5(cos#-cos# ) 6(^2 

P2q 



-Me- 



(37) 



xpI 
2m: 



pI 



X 2 



1 + 



xMlc 

q 2 , 



q 2 



P\ 



P 2 ± 



2M 2 * 1 2M 2 * ' 

where the new variables appearing above are defined by the relations 



M* 

X = 1 ~ AT ' c = <7o + ^2-t/4- 



2 Ml 



(38) 



(39) 



The factor U 2 — U4 is the potential energy gained in converting a particle of species "2" to 
a particle of species "4". When the intial and final state baryon masses are different, the p 2 
integration is restricted to the interval p\ > p\ > p 2 _. Note also that the effective chemical 
potential \i 2 — U 2 replaces /12. With these changes, the response function incorporates the 
effects of strong interactions at the mean field level. We can now generalize the definition 
of which appears in the noninteracting response function, with £±: 



In 



1 + exp((e± - n 2 + U 2 )/T) 



1 + exp((e± + g - + U 2 )/T)_ 
Collecting these modifications together, and using the relation in Eq. 



we have 



S(qo,q) 



M*MIT 
nq 



exp(— z) 



(40) 



(41) 
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For the charged current, modifications due to interactions are twofold. First, the difference 
in the neutron and proton single particle potentials appears in the response function and 
also in fi. Second, the response depends upon the nucleon effective mass. Since both and 
M* are strongly density dependent, the opacities are significantly altered from those for the 
noninteracting case at high density. 

Finally, note that we can make the nondegenerate, elastic, and degenerate approxima- 
tions in analogy to those found for the noninteracting case. The expression for the degenerate 
approximation is nearly indentical to that for the noninteracting case, because the condition 
e - > A*2 — U 2 is equivalent to q > \p F2 — p Fi \. We must only replace M 2 in Eq. (GBf) by 



M2M4. The expression for the nondegenerate approximation is identical to that for the 
noninteracting case (Eq. (|31~1) ), because for nondegenerate nucleons, the number density is 
proportional to (MjfT) 3 / 2 exp((/i2 — U 2 )/T). Finally, in the elastic case, we need only replace 
/i 4 — /i 2 in Eq. (|4j) with /i 4 — /i 2 + U 2 — U^. 

For the neutral currents, since the initial and final state particle labels are identical for 
both leptons and baryons, we have the simplifications 



2 \2 



z=Tf> ^ = ^ e - = 2q^{ q °-2M*) ' e + = °°- (42) 
Thus, = — z, and one finds a result formally similar to the noninteracting expression: 



S{qo,q) = 

nq 



— 1 + 

exp(— z) V 



(43) 



Fig. 8 shows the composition and the chemical potentials versus the density for charge- 
neutral matter containing interacting nucleons, electrons, and trapped neutrinos in beta 
equilibrium. While these results are qualitatively similar to the case of noninteracting nu- 
cleons, interactions lead to lower values of Y Ue and [i„ e , and larger values of jl, especially at 
high densities. The influence of these changes on the neutrino mean free paths are depicted 
in Fig. 9, using E\ = fi Ue as before for the trapped-neutrino case, and the dynamic form 
factor Eq. fl4~l~|). Shown for comparison are the three limiting behaviors for the degenerate, 
nondegenerate, and elastic approximations discussed in connection with the noninteracting 
case. 

The counter-intuitive result, that the mean free paths increase at high density, is chiefly 
due to the use of nonrelativistic kinematics coupled with the behavior of the Landau effec- 
tive mass of the nucleons, which are displayed in Fig. 10. The cross section is proportional 
to M2M4, and the dropping effective masses more than compensate for the higher neutrino 
chemical potential at high densities. However, the net effect of a decreasing cross section 
at high density may be an anomalous result of using a nonrelativistic theory in a situation 
in which the nucleons are at least partially relativistic. This also explains why the elastic 
approximation, which is applicable in the case when the mass of the absorber is large, pro- 
vides a much poorer representation of the exact results than in the noninteracting situation. 
Related to this is the fact that the high density behavior of the nonrelativistic Skyrme-like 
EOS becomes dominated by the w 8 / 3 density-dependence in the energy density. This leads to 
eventual acausality, i.e., the speed of sound exceeds the speed of light. While it is possible to 
avoid this behavior in nonrelativistic potential model approaches by including more general 
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forms of momentum dependence (see for example, Ref. |48|]), a nonrelativistic description of 
matter at very high densities is difficult to justify at a fundamental level. In addition, the 
calculation of the opacity is then greatly complicated by the fact that the energy spectrum 
no longer resembles that of a free gas. In the next section, we consider field-theoretical 
models with relativistic kinematics which do not have the above problems. 

In the neutrino-free case, interactions change the composition significantly. The proton 
fraction is larger when compared with the neutrino-free noninteracting case. Quantitatively, 
the proton fraction is sensitive to the density dependence of the symmetry energy, which for 
the Skyrme interaction we chose is nearly linearly increasing with density. The composition 
and chemical potentials are shown in Fig. 11. We shall first focus on the charged current 
reaction. At zero temperature, the direct Urea reaction v + n — > e~ + p is kinematically 
possible for low energy neutrinos at and above a threshold density when the proton fraction 



exceeds 1/9 |45 |. In the presence of muons, the proton fraction at threshold density is 
slightly larger. The threshold density for the model chosen here occurs at n = 1.5n (for a 
free gas, by comparison, the threshold density is n = 73n !). 

In Fig. 12 the neutrino mean free paths are calculated for the neutrino-free case. At low 
temperature, the absorption mean free path increases sharply below the threshold density 
(1.5no). The mean free paths show similar behavior to that of the neutrino-trapped case 
shown in Fig. 9, although in general they are about 3 times smaller. This factor can largely 
be understood by examining the last three terms in the noninteracting expression Eq. ( ^6|) . 
In the neutrino-trapped case in which Ei « fi u these terms yield [i e (YL = 0.4)7r 2 /2. In the 
neutrino-free case in which ji v = and E\ = 3T, these terms yield \i e {Y v = 0)(n 2 + 9). The 
ratio of these two cases is thus approximately fi e (Yi = 0.4)/4.5/i e (Y"^ = 0) ~ 0.3 using values 
for [i e from Figs. 8 and 11. 



V. RELATIVISTIC INTERACTING BARYONS 

At several times the nuclear equilibrium density uq, the Fermi momentum and effective 
nucleon mass are both expected to be comparable. Thus, a relativistic description may be 
more appropriate. Relativistic local quantum field-theoretical models of finite nuclei and 



infinite nuclear matter have had some success ||35|| , albeit with rather more schematic inter- 
actions and with less sophisticated approximations than their nonrelativistic counterparts. 
Appendix B contains a description of this approach in which the set of baryons B has been 
augmented to include strangeness-bearing hyperons. 

The usual starting point for relativistic field theory calculations is the mean-field approx- 
imation. In this approximation, the presence of baryons generates nonzero average values 
of the meson fields. Baryons move independently in these self-consistently generated mean 
fields. In a model of baryons Yukawa coupled to vector (^), isovector (b^) and scalar (a) 
mesons, the latter simulating correlated pion-exchange, the single particle energies are 



Eb = VP 2 + M$ + guBWo + g P Bh B b = VP 2 + Mg + U, (44) 

where M B = M B —g aB a are the Dirac effective masses, a, uj , and b are the average values of 
the meson fields, and g aB , g^B-, and g pB are the strong interaction couplings of the different 
meson fields to baryons, and t^s is the third component of isospin for the baryons. This 
equation defines the effective potential U. For more details, see Appendix B. 
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We turn now to calculate the cross section for absorption of electron neutrinos in a 
multi-component system described by a Walecka-type effective field-theoretical model at 
the mean field level. In this model, the nucleons become increasingly relativistic due to 
the rapidly dropping nucleon effective masses. The nonrelativistic approximation for the 
hadronic current discussed thus far cannot be justified, since, at high density, terms of 
order p/M* are of order unity and must be retained. It is convenient to express the angles 
appearing in Eq. (0) in terms of energy and momentum transfers (qo,q), and to define 
the response functions in terms of appropriate current-current correlation or polarization 
functions P7| . To see how this is accomplished, we start with a general expression for the 
differential cross section [ [49| , |50[ 



d 3 a G\ E 3 r ' 1 ' ' ' 



V d?VL 3 dE 3 128tt 2 Ex 



exp I 



f 3 (E 3 )) Im (L^H*) , (45) 



where the incoming neutrino energy is E\ and the outgoing electron energy is E 3 . The 
factor [1 — exp((— q — H2 + Ha) /T)}~ 1 arises due to the fluctuation-dissipation theorem, since 
particles labeled '2' and '4' are in thermal equilibrium at temperature T and in chemical 
equilibrium with chemical potentials fi2 and /14, respectively. The final state blocking of the 
outgoing lepton is accounted for by the Pauli blocking factor (1 — f 3 (E 3 )). The lepton tensor 
L af3 is given by 

L a/3 = 8[2k a kP + (k ■ q)g aP - (JfeV + q a k P ) T ie^Wq"] (46) 
The target particle retarded polarization tensor is 



Ml^teri. " T ^-"" lmn tf , (47) 



where H a p is the time ordered or causal polarization and is given by 

= -i J ^Tr [T(G 2 (p)J a G 4 (p + q)J p )] . (48) 

Above, is the incoming neutrino four-momentum and g M is the four-momentum transfer. 
In writting the lepton tensor, we have neglected the electron mass term, since typical electron 
energies are of the order of a few hundred MeV. The Greens' functions G{(p) (the index i 
labels particle species) describe the propagation of baryons at finite density and temperature. 
The current operator J M is 7^ for the vector current and 7^75 for the axial current. Given 
the structure of the particle currents, we have 

rW = V 2 nr, + A 2 U% - 2VAU V J • (49) 

For the vector polarization, {J a , Jp} {la,lp}, for the axial polarization, {J a ,Jp} 
{la 75, 7/375} and for the mixed part, {J a ,Jp} :: {7a7s,7/3}- Using vector current con- 
servation and translational invariance, nig may be written in terms of two independent 
components. In a frame where g M = (q , \q\, 0, 0), we have 



a 2 

IT T = Ii v 22 and n L = -f^Ii v 0Q 

m 
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The axial current-current correlation function can be written as a vector piece plus a cor- 
rection term: 



n 



(50) 



The mixed, axial current-vector current correlation function is 



n 



VA 



ie 



VA 



(51) 



The above mean field or Hartree polarizations, which characterize the medium response to 



the neutrino, have been explicitly evaluated in previous work 0,0] for the case of neutrino 
scattering. Since the structure of the charged current is similar to that of the neutral current, 
we may write 



d 3 



a 



V d?VtdE v 



Cl #<fi(l-/a(^)) 



A 



16tt 3 E l 
4:EiE 3 + q\ 
2q 2 



exp 



"<?0 - [fJ-2 - A*4) 



n -1 



T 



B = E 1 + E 3 



[ARt + R 2 + BR 3 ] , 



(52) 



From Eq. (^), we see that the response of a relativistic system to the charged or neutral 
current probe may be written in terms of three response functions R±, R 2 , and -R3. (In 
contrast, the response in a nonrelativistic system, Eq. (|14D , is characterized by a single 
isospin density response function. This is due to the fact that when the baryon velocity 
Vi ~ 1, the angular terms in Eq. (0) cannot be dropped.). The various response functions 
required have been studied in earlier works for neutral current reactions [P7,p3]. Their 
generalization to the case of charged current reactions are: 



R x = (V 2 + A 2 ) [Im nf (g , q) + Im n£(g , q)} 
R 2 = (V 2 + A 2 ) Im n£(g 0) q) - A 2 Im n^(g , q) 
R 3 = 2VA Im n^ A (g ,g). 



(53) 
(54) 
(55) 



These response functions have been written in terms of the imaginary part of the polarization 



functions, whose causal components are given in Ref. f51"| , |52|| for symmetric nuclear matter. 
We present extensions of these results to asymmetric matter, and, in particular, to unlike 
particle- hole excitations. For space like excitations, q 2 < 0, they are given by 



i m n i(M>2W ^ e; 2 -\ p \ 2 C qs 2 8 

1 (2tt)3 E;E* p+q 
d 3 p q 2 J2 - |p| 2 (l - cos 2 



Im n T (g , £)= fi- 



lm n A (go, q) = 2vr J 
Im U VA (q ,q) = 2tt / 



(2tt) 3 
d 3 p 



(2vr)3 E;E. t 
d 3 p 



m( e 



e 



p^p+q 



(2tt)3 \q*\E;E* p+q 



e 



(56) 
(57) 
(58) 
(59) 



In the above, 



17 



e = f(e;, E; +q )[6( go - (E p+q - e p )) + 6( qo - (e p - E p+q )] 
F(E;,E; +q ) = f 2 (E* p )(i-fr(E; +q )) 

e p = e: + u 



e; = ^/\p\* + m?, ^ ^, 

The particle distribution functions fi(E) are given by the Fermi-Dirac distribution functions 

1 



(60) 
(61) 

(62) 



me; 



1 + exp((E* p - Vi)/kT) ' 
where v is the effective chemical potential defined by 

Vi = fa - Ui = ^ - (guBiUo + t 3Bi g P Bibo) 



(63) 



(64) 

The single particle spectrum in Eq. (f44|) is discussed in Appendix B. The angular integrals 
are performed by exploiting the delta functions. The three dimensional integrals can be 
reduced to the following one dimensional integrals: 

,2 



Im n L (g , q) = jf dE [(E + g /2) 2 - |g| 2 /4] 

x [F(E,E + q ) + F(E + q ,E)} 
Im n T (g , q) = J^L ^ dE [(JET + g /2) 2 + |g| 2 /4 + |9|X7#)1 



Air\q 

x [F(E,E + q ) + F(E + q ,E)} 



Im U A {qo,q) 



Im U VA (q ,q) 



M* 



2Ti\q\ 



57T|g 



roo 

^ dE [F(E,E + q ) + F(E + q ,E)} 

roo 

- J dE [2E + g ] [F(E, E + q ) + F(E + q , E)] . 



The lower cut-off e_ arises due to kinematical restrictions and is given by 



-8^ + 1 



2 



where 



qo = qo + U 2 -U i , 6=1 



MX 2 
[ 

(f ~ Qo 



Mf - u.r 

q 2 -ql 



It is convenient to reexpress the polarization functions as follows: 



Im n£(g ,g) 
Im ITf (g , g) 

Im n«(g ,g) 
Im n^ A (g ,g) 



h + qoh + 



h + qoh+{^ + q j + Mf^I 
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(65) 

(66) 
(67) 

(68) 

(69) 

(70) 

(71) 
(72) 

(73) 
(74) 
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where we used the one- dimensional integrals 



tanh 



go + O2 - At 4 ) v 
2T 



^ £ n [F(E, E + q )+ F(E + g , E)] . 



(75) 



These integrals may be explicitly expressed in terms of the Polylogarithmic functions 

r* Lin-Ax) 



LiJz) 



dx . 



x 



Li 1 (x) = ln(l — x) 



(76) 



which are defined to conform to the definitions of Lewin |53| . This Polylogarithm represen- 
tation is particularly useful and compact: 



h = Tz\\ 
h=T 2 z 
h=T 3 z 



T 

(M2 - W 
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2 z zT y 

- z^ + — + 2— + 2 
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+ 



\ T 2 T 3 3 z Tz T 2 z J 

where z = (go + (/i 2 — fii))/T and the factors £ n are given by 

£ n = Li n (-ai) - Li n {-oi 2 ) , 

with 

oti = exp ((e_ - /i 2 + t^/T) , a 2 = exp ((e_ + g - /i 4 + t/ 2 )/T) 



We note that the nonrelativistic structure function for neutral current scattering, Eq. (|43| 
is, aside from the factor MJ T/irq, equal to Iq since £1 = 

The total cross section is the double integral in (g , g) space: 



(77) 
(78) 
(79) 

(80) 
(81) 
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(82) 



Eq. (^) allows us to calculate the cross section per unit volume, or equivalently the inverse 
mean free path, consistently with the relativistic field-theoretical model in the mean field 
approximation. This naturally incorporates the effects of strong interactions, Pauli blocking 
of final state particles and the contribution of relativistic terms to the baryon currents. In 
the case of neutral currents, some of the terms in the above are slightly simplified: 



go 

T 



A*2 — A*4 



% q 
' 2 2\ 



Mf 

' ql 



(83) 



We note that other expressions for the cross sections, based on a similar formalism, 
have been derived earlier in the literature [^,^]. However, negative values of go were not 
taken into account in Ref. |27j and the response functions used in Refs. [p7||33f inadvertantly 
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omitted the factor (1 — exp(— z)) in the denominator of Eqs. (|52|). While the qualitative 
results in Refs. p7|j33[| remain unchanged, they are quantitatively affected. 

In Fig. 13, the individual particle fractions and the relevant chemical potentials are shown 
for a lepton fraction Yi = 0.4. As in the case of the potential model, interactions in this 
model lead to larger proton fractions compared to the noninteracting case. These results are 
almost indistinguishable from those of the potential model. The absorption (upper panels) 
and scattering (lower panels) mean free paths, calculated using Eq. QS2]), are shown in Fig. 14 
for three different temperatures, T = 5, 30, and 60 MeV, as a function of baryon density. At 
low temperatures, the neutrino mean free path has relatively little variation with density due 
to the dropping baryon effective masses at high density. Nevertheless, the mean free paths 
eventually decrease with increasing density, unlike in the potential model. We attribute this 
to the inclusion of relativistic kinematics. 

The results for the neutrino-free case are shown in Figs. 15 and 16. Qualitatively, the 
results have similar density and temperature behaviour as in both the neutrino-trapped 
case and in the nonrelativistic case. The upper panels show the absorption mean free paths 
and the lower panels show the scattering mean free path. The threshold density for the 
charged current reaction to be kinematically allowed at zero temperature in this model 
occurs at 1.7n . This accounts for the sudden increase in absorption mean free path below 
this density at low temperatures. At higher temperatures, this threshold-like behaviour is 
much less pronounced, neutrino absorption mean free path dominate the total opacity for 
electron neutrinos. From the leftmost panels in Fig. 16, we see that even at T = 5 MeV, 
despite kinematical restrictions at low density the absorption reaction always dominates over 
the scattering reaction. 

VI. COMPARISON WITH OTHER WORKS 

Calculations of protoneutron star evolution have been performed by several groups § . 
At supra-nuclear densities, these groups used for the most part neutrino cross sections as 
originally described in Burrows and Lattimer Q. Those cross sections were generated by 
interpolation among limiting formulae for degenerate and nondegenerate matter and neutri- 
nos analogous to the limiting formulae given in Sec. Ill for noninteracting, nonrelativistic 
baryonic matter. A direct comparision of those results with the results obtained here for 
arbitrary degeneracy and including interactions is therefore not straightforward. First, the 
limiting expressions should include the effects of interactions. Second, the results will be 
sensitive to the interpolation algorithm, especially the interpolation parameter. 

In addition, relativistic effects become increasingly important with density in a field- 
theoretical description, since the nucleon effective mass decreases with increasing density. 
Quantitatively, the effects arising due to the relativistic structure of the baryon currents is 
small and is typically of order Eu/M* . However, relativistic kinematics introduces important 
corrections (of order pf/M*) in the degenerate limit. To facilitate an illustrative comparision 
of our results with other results often used in the literature, we study the neutrino mean 
free path in the trapped regime. In this case, the dominant neutrino cross sections are those 
due to absorption (on neutrons) and scattering off neutrons and protons. The appropriate 
limiting expressions for the absorption and scattering cross sections for the degenerate case 
become 
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where 



Ppj + M* 2 replaces M 2 appearing in the formulae for nonrelativistic, nonin- 
teracting nucleons Eq. ( p6|) (z = n,p). For scattering off neutrons and protons, cy and ca 



are appropriately chosen. The nondegenerate cross sections n o^td j are unchanged 
from Eq. (0). 

We compare the results for the total cross section obtained by interpolation as opposed 
to the exact integrations (<J rp i) for the case of neutrino-trapped matter with Y L = 0.4. For 
this exercise, we chose a field-theoretical model and therefore used the limiting formulas 
for the cr^'s and ct^d's just described. The total cross section is the sum of those from 
absorption on neutrons and scattering off neutrons and protons. As an illustration, we first 
employ the particular interpolation scheme of Keil and Janka M : 
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where X{ = max[0,^j] where 77, is the degeneracy parameter. For a relativistic model, the 
appropriate degeneracy parameter is rji = {i>i — M*)/T which can be seen by reference to 
Eq. (|63D . Fig. 17 (lower left panel) shows the ratio of our exact integration compared with 
o\ as a function of density and temperature and assuming E v = \i v . This interpolation gives 
reasonable cross sections in the nondegenerate limit but is poor in the degenerate limit. 
This behavior is easily understood by examining the ratio cr^ D /o"^, which for E u = n u is 
proportional to nE 2 / \E*p n T 2 /i e ) . We can neglect the relatively weak density dependence of 
Ep , since the drop in M* is compensated somewhat by the increase of kp n , and note that in 
the neutrino-trapped case both E v = [i v and /z e scale as n 1 ^ 3 . Thus, ct^d/pd ~ n^^/T 2 ~ 
(Epn/T) 2 = rj 2 where Ep n is the nonrelativistic neutron Fermi energy. In the degenerate 
limit, therefore, the <x^ D contribution can still dominate that from a^. This may also be 
inferred by comparing with the ratio <J rv ij<JND i n the degenerate regions as shown in the 
upper left panel of Fig. 17, where ctjvd — °nd + °1vd- 

These results suggest that a better interpolation formula is 
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The ratio of the exact to these interpolated results are shown Fig. 17. The interpolation now 
correctly goes to the required limits; for reference the exact result is compared directly with 
the two limiting forms in the upper left (nondegenerate) and upper right (degenerate) panels. 
Nevertheless, significant errors for intermediate degeneracies still exist, which underscores 
the importance of the relatively simple exact formulae we have found. There is no significant 
need to use interpolations any longer. 
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VII. MULTI-COMPONENT ENVIRONMENT 



Strangeness-bearing components can appear in dense matter either in the form of hy- 
perons, kaon condensate, or quarks. In this section, we investigate the effects of multi- 
components on the neutrino mean free paths by concentrating on the possible presence of 
hyperons. We present specific results for a relativistic field-theoretical model in which the 
baryons, B, interact via the exchange of cr,u, and p mesons [p5[| . The relevant details are 
presented in Appendix B. 

For a typical case, the particle fractions in matter containing hyperons are shown in the 
upper panels of Fig. 18 for different temperatures in the neutrino-trapped case, Y L = 0.4. 
The lower panels show the behaviour of the relevant chemical potentials. Note that the 
neutrino chemical potential increases rapidly with the onset of hyperons in matter (compare 
with the results in Fig. 12) due to the increase in the positive charge in the system. This is 
typical for matter in which strangeness, due to any source, appears. Consequently, in matter 
with hyperons neutrino energies are somewhat larger than those in nucleons-only matter. 

The absorption mean free paths in matter containing hyperons are shown in Fig. 19. 
The various reactions that contribute to the total absorption opacity are calculated using 
Eq. ( jS2j) by using the appropriate chemical potentials and neutrino coupling constants. 
Strangeness-changing reactions are Cabibbo suppressed and hence their contribution to the 
total opacity is negligible. The relative importance of the various reactions are shown in 
the upper panels. Pauli blocking and kinematic restrictions account for the threshold-like 
structure seen at low temperatures, particularly for the reaction u e + S~ — ► A + e~ . This is 
analogous to the threshold behavior seen for the reaction v e + n —>■ e~ + p in hyperon-free 
matter. The kinematical restrictions on reactions involving baryons require that the initial 
and final state baryon Fermi momenta not be vastly different. When too large a difference 
exists, the phase space is highly suppressed. At higher temperatures, these kinematical 
restrictions are relaxed. The higher neutrino energy and the lower degeneracy of the baryons 
account for the qualitative trends of the results at the higher temperatures T = 30 MeV 
and T = 60 MeV. The solid lines in the lower panels show the net mean free path from all 
contributions from absorption reactions. For comparison, the net mean free path in a model 
without hyperons is shown by the dashed curves. In general, the presence of hyperons has 
the effect of decreasing the neutrino mean free path. 

In Fig. 20, the scattering mean free paths are shown. The upper panels show the indi- 
vidual contributions to the total scattering mean free path due to the various reactions of 
interest. In the lower panels, the net mean free paths in hyperonic matter (solid curves) 
are compared with those in nucleons-only matter (dashed curves). The appearance of hy- 
perons again decreases the neutrino mean free path. The overall scattering mean free path 
is less than 50% larger than the absorption mean free path, so that scattering provides an 
important contribution. 

In summary, the decrease in mean free paths in a mult i- component environment may 
be understood by noting that: (1) A larger neutrino chemical potential results in larger 
neutrino energies, (2) hyperons decrease the degeneracy of baryons, and (3) hyperons provide 
additional channels of scattering reactions to occur. 

When neutrinos carry no net lepton number, they do not play a role in determining the 
composition of matter. In neutrino-free matter, the hyperons appear in larger numbers and 
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at threshold densities which are lower than in neutrino-rich matter. In Fig. 21, particle 
fractions (upper panels) and the electron chemical potentials (lower panels) in neutrino-free 
matter with hyperons are shown. Since the neutrinos in this case are thermal, we present 
results for absorption and scattering mean free paths calculated for E v = 3T. 

In Fig. 22 are shown the relative contributions of the dominant charged current reactions 
(upper panels) and the total neutrino absorption mean free paths (lower panels). At low 
temperatures, the reactions show the expected threshold-like behavior. The total absorption 
opacity due to all possible reactions (solid curves) is shown in the lower panels. For com- 
parison, results in matter without hyperons is also shown by the dashed curves. As in the 
neutrinos-trapped case, the neutrino mean free paths are smaller in matter with hyperons, 
albeit by a smaller amount. Scattering off hyperons, shown in Fig. 23, clearly leads to 
important modifications to the neutrino mean free path in the neutrino-free case. 

VIII. NEUTRINO TRANSPORT IN AN EVOLVING PROTONEUTRON STAR 

A. Introduction 

One of the most important applications of dense matter opacity calculations is the en- 
vironment of a newly formed neutron star. A protoneutron star is formed subsequent to 
the core bounce of a massive star in a gravitational collapse supernova. Its early evolution 
has been investigated in Refs. |§-§|. Detailed studies of the dynamics of core collapse and 
supernova indicate that within milliseconds of the shock wave formation the formerly col- 
lapsing stellar core settles into nearly hydrostatic equilibrium, with a relatively low entropy 
and large lepton content. 

The entropy per baryon, s, is about one or less (measured in units of Boltzmann's 
constant), which corresponds to a temperature of about 10-20 MeV. The electron lepton 
fraction Yl = Y e + Y Ve at bounce in the interior is estimated to be about 0.4. The z/ e 's 
formed and trapped in the core during collapse are degenerate with a chemical potential 
of about 300 MeV. In addition, because no fi- or r- leptons were present when neutrino 
trapping occurred, the net numbers of either fi or r leptons is zero: e.g., Y v = —Y^ « 0. 
The /i and chemical potentials in beta equilibrium are related by fi^ — fj,^ = fi e — fi Ue , each 
side of which has a value of order 100 MeV. Unless /i M > m M c 2 , however, the net number of 
/i's or z/^'s present is zero. 

The temperature and entropy increase beyond about 0.5 M Q from the center because of 
shock heating H. During the early deleptonization or neutrino loss phase, the stellar interior 
gains entropy because of resistive neutrino diffusion, and the regions near the protoneutron 
star surface lose entropy because of neutrino losses. Eventually, this reverses the positive 
temperature and entropy gradients in the interior. On time scales of about 10-15 s, the 
center reaches a maximum entropy of about 2 and a maximum temperature of 40-60 MeV. 
This time coincides with the loss of virtually all the net trapped neutrino fraction Y Ue in the 
interior, although there are still considerable numbers of neutrino pairs of all flavors present 
in thermal equilibrium. This time, therefore, marks both the end of the deleptonization 
phase and the onset of the cooling phase of the protoneutron star. 

The baryonic composition of dense matter is greatly affected by the degree to which 
neutrinos are trapped. Thus there exists an unambiguous compositional difference between 
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the initial deleptonization (neutrino-trapped) and cooling (neutrino- free) phases |47[]. This 



difference could be enhanced if strangeness, in the form of hyperons, kaons or quark matter, 
is considered. Neutrino trapping delays the appearance of the strange matter components 
to higher baryon densities. This implies that during the early deleptonization phase, matter 
may consist mostly of nonstrange baryons, except possibly at the very center of the star. 
But the cooling phase may be characterized by the presence of a substantial amount of 
strange matter, due to both the decreasing threshold density for their appearance and to 
the increasing central density of the star. It is therefore of considerable interest to examine 
the behavior of neutrino opacities in matter characterized by these different thermodynamic 
conditions. 



B. Semi-analytic treatment of neutrino transport 



The neutrino transport in the bulk of the interior may be treated in the diffusion ap- 
proximation, since neutrinos are very nearly in thermal equilibrium due to the large weak- 
interaction rates and the typical mean free path is very small compared to the stellar radius 
until after about a minute when the mean neutrino energy becomes very small. In the diffu- 
sion approximation, neglecting the effects of general relativity, the rate of change of electron 
lepton number is related to the electron neutrino number gradient by 
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where n is the baryon number density, 
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are the electron neutrino and antineutrino number densities, respectively, at energy E v , 
f u = (1 + e( E "~ ,J "'" T )~ 1 and f p = (1 + e^ Eu+tJ,u >' T )~ l are the neutrino and anti-neutrino 
distribution functions, and Yl = Y v + Y e is the total number of leptons per baryon. The net 
neutrino mean free path is due to both absorption and scattering: 
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where o a and a s are the absorption and scattering cross sections, respectively, and the factor 
1 — fv{E v ) accounts for the inverse process |T7[. A similar relation exists for the antineutrino 
mean free path, but uses f„ and the appropriate absorption and scattering cross sections 
{o'ai &s) for u e s. Eq. (RBI) can be rewritten in terms of the diffusion coefficients 
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The diffusion constants are related to the conventional Rosseland mean free paths by 
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The diffusion constants D2, -D3, and D% only contain contributions from electron-type 
neutrinos, since only these are associated with changes in electron-neutrino number. 

During the deleptonization phase, the neutrinos are mostly degenerate. Thus, the diffu- 
sion of antineutrinos can be essentially ignored and the neutrino properties are essentially 
only a function of \x v . For nearly degenerate matter in beta equilibrium, one can show that 
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where the subscript o indicates values at the beginning of deleptonization. For example, for 
fi U)0 ~ 200 MeV, (dY L /dY u ) Q ~ 3. It is instructive to explore the transport behavior for 
various assumptions concerning the energy dependence of the mean free path. Supposing 
that \ V (E V ) = \ (n,T)(fi Ut0 / E u ) m , one finds in the degenerate limit 
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Therefore, assuming that during the deleptonization phase the remnant is nearly hydrostatic, 
Eq. (|92"D can be written in terms of the neutrino chemical potential as 
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Note that the terms involving the temperature gradient in Eq. ( |92f) vanish in the degenerate 
limit irrespective of the value of m. The temperature now only enters through A G (n, T). 
Inasmuch as the temperature in the interior of the neutron star varies with time only by 
a factor of about two during the deleptonization, while //„ varies by a much larger factor, 
a simple understanding of deleptonization can be obtained by treating A Q as a constant in 
both space and time. 

As shown in Ref . j|7| , approximate solutions of Eq. (^Bp can be found by separating the 
time and radial dependences in \i v : 



fJ> v = fj. v , ij>(t)il}(r) 



(97) 



with 0(0) = 1 and i/)(0) = 1. This separation is justifiable since the remnant is essentially 
hydrostatic. The functional forms of (f) and ip, and the eigenvalue of the solution, depend 
upon the value of m. In Ref. [|7|], it was assumed that m = 2, in which case <j) decreases 
linearly with time: <fi — 1 — t/r. If, instead, m — 0, then <ft varies like (1 + t/r)- 1 . In the 
case of m = 1, varies like exp— (t/r). In each case, r is a deleptonization time and is 
proportional to R 2 /c\ where R is the stellar radius. In detail, one finds 
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|54| . For the cases 
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where is the Lane-Emden radial eigenvalue for index n = 
m = 0, 1 and 2, one finds n = 2/3, 1, and 2, and £„ i = 2.871, 
low, we examine the behavior of D 2 with respect to //„, for fixed temperature, in the interior 
of a star at the onset of deleptonization, as it appears that the details of deleptonization 
fundamentally depend upon this behavior. 

As the net electron-neutrino fraction disappears, the number diffusion equation becomes 
irrelevant and one should consider the energy diffusion equation [07 
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Here, s is the entropy per baryon and e v j = n v /E v are the energy densities, at energy E u , 
of neutrinos (or antineutrinos) of species I. Similarly, X u ^ represents the net mean free path 
of each neutrino and antineutrino species. This equation assumes the matter to be in beta 
equilibrium. The energy densities and mean free paths, in contrast to the deleptonization 
situation, contain important contributions from all three types of neutrinos. The A's of 
electron-type neutrinos are due to both absorption and scattering, but those of tau-type 
neutrons have only scattering contributions. Those pertaining to muon-type neutrinos, 
besides having scattering contributions, may also have absorption contributions if muons 
are present, which is generally the situation only when Y v < 0.02 above nuclear densities. 
In the neutrino-free case, setting \i v = in Eq. fl99|) results in 
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As is the case with the mean free paths, the diffusion coefficients D±t have contributions 
from all three types of neutrinos. Since the baryon matter, which dominates the specific 
heat, remains fairly degenerate throughout the cooling epoch, the entropy varies practically 
linearly with the temperature. In the nondegenerate approximation for the baryons and the 
neutrinos, the relativistic cross sections are expected to vary as E^n, so that oc T 3 /n. 
However, the baryons are degenerate so we should expect the cross sections to be modified 
by an extra degeneracy factor of T j \i n . Therefore, we can anticipate that the diffusion 
constants D^e will behave as T 2 \i n jn. A separation of Eq. (|100|) into temporal and spatial 
variations reveals that the temperature in the star is then expected to be roughly linearly 
decreasing with time. The predicted temperature profile will depend somewhat upon the 
density profile of the star; assuming the density to be approximately constant, one finds that 
the temperature profile is that of an n = 2 Lane-Emden polytrope. This is different from the 
predictions in Ref . [0 , which, however, did not account for the degenerate nature of baryonic 
matter during the cooling stage. However, numerical simulations have demonstrated the 
qualitative correctness of these remarks. 



C. Results 

Fig. 24 shows D 2 plotted versus density for beta equilibrium matter with Yl = 0.4 for 
temperatures between 5 and 30 MeV for the field-theoretical model adopted earlier. Both 
nucleons-only matter and matter in which hyperons are assumed to appear are shown. One 
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sees that both nucleons-only matter and matter containing hyperons have relatively small 
variations of D 2 with density at fixed temperature for n > 1.5n s . In the limit of degenerate 
matter and neutrinos, this results from the asymptotic behavior of the cross sections at the 
Fermi surface, which is a oc T 2 /z e /z n /i p in the relativistic case. (In the nonrelativistic case, 
see Eqs. (pip and (|27j).) Thus, D 2 is expected to behave like /i^/(T/i e /x n /x p ) in this case, and 
this has relatively little density dependence at high densities since the increase in \i e \i n \x v is 
compensated by the increase in This implies that the case m = 2 could approximate the 
deleptonization stage of a protoneutron star, which implies that [i v should decrease linearly 
with time. Since D 2 is roughly proportional to 1/T, one has that A G is roughly proportional 
to 1/T 2 . The rising temperature during deleptonization will modify the linear decrease of 
\x v . However, numerical simulations show that the linear behavior is approximately correct. 

The effect of hyperons appearing during the deleptonization phase at first does not appear 
to be appreciable, especially if hyperons appear only relatively late in the evolution. In the 
presence of hyperons, the diffusion constant D 2 decreases from its value for the nucleons-only 
case for a given temperature. However, the appearance of hyperons is accompanied by a 
general temperature decrease because of the addional degrees of freedom which increase the 
heat capacity of the system. This results in a very small net change for deleptonization times. 
More important effects from hyperons can be expected at the late stages of deleptonization 
and during the cooling epoch following deleptonization. The total -D 4 = J2e D^g is shown as 
a function of density for several values of the temperature between 5 and 30 MeV in Fig. 25. 
As suggested in the previous section, an overall T 2 behavior is observed. The large (factor 
of 2) decrease in in hyp eron-b earing matter relative to nucleons-only matter for a given 
T would, at first glance seems to imply longer diffusion and cooling times. However, the fact 
that the temperatures are smaller in the presence of hyperons somewhat compensates for 
this effect and preliminary calculations show that there is an overall increase in cooling 
times for models which contain hyperons. Thus, there are important feedbacks operating 
between the stellar structure, the EOS and the opacities that may only be addressed in 
detailed simulations. 



IX. SUMMARY AND OUTLOOK 

In summary, we have calculated both charged and neutral current neutrino cross sections 
in dense matter at supra-nuclear densities. We have identified new sources of neutrino 
opacities involving strange particles and have computed their weak interaction couplings. 
The weak interaction cross sections are greatly affected by the composition of matter which 
is chiefly determined by the strong interactions between the baryons. We have, therefore, 
performed baseline calculations by considering the effects strong interactions on the in- 
medium single particle spectra, which also determines the composition through the EOS. 
The formalism we have developed allows us to calculate the cross sections efficiently for 
matter at arbitrary matter degeneracy. From our results, various limiting forms used earlier 
in the literature are also easily derived. This, in addition to providing valuable physical 
insights, allows us to assess the extent to which the various approximations are valid in both 
free and interacting matter. 

To explore the influence of baryonic interactions on the neutrino cross sections, we have 
examined both nonrelativistic potential and relativistic field-theoretical models that are 
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commonly used in the calculation of the EOS. We sought to identify the common features 
shared by these models. At the mean-field level in both relatively simple structure 

for the single particle spectrum leads to analytic expressions for the response function, which 
in turn facilitates a first study of the role of strong interactions on the neutrino cross sections 
at temperatures of relevance in astrophysical applications. Specifically, when the effects of 
momentum dependent potential interactions can be treated adequately in the effective mass 
approximation (i.e., only the quadratic term in momentum is retained), analytic expressions 
obtained for noninteracting matter are straighforwardly extended to include interactions. 
Effective field-theoretical models at the mean-field level offer a similar opportunity since 
the single particle specturm is that of a free Dirac spectrum, but with a density dependent 
effective mass. Note that in both cases, energy shifts in the spectrum arising through density 
dependent interactions are naturally included. Investigating the effects of more complicated 
momentum dependent interactions will be necessarily more involved. We have also calculated 
neutrino cross sections in matter containing hyperons and assessed their influence on the 
total opacities. 

We have examined the role played by neutrino opacities in determining the deleptoniza- 
tion and cooling times of a newly born neutron star as it evolves in time. We showed 
analytically how these times are related to the relevant diffusion constants and, thus, the 
opacities. Although the main physical issues involved are clarified in such an approach, a 
quantitative assesment must await calculations using a more complete protoneutron star 
evolution code. The results of such a calculation will be reported elsewhere []55 . 

Some important aspects of our work are: 

1. An exact and efficient calculation of the phase space valid for arbitrary degeneracy of 
asymmetric matter for both nonrelativistic and relativistic interacting baryons. The 
simplicity of the resulting formulae may be especially useful in the calculation of elec- 
tron scattering, for example. 

2. With the formalism presented here, it is no longer necessary or desirable to em- 
ploy "generic" (i.e., EOS independent) neutrino opacities in astrophysical simulations. 
Rather, tables of opacities can be constructed directly from quantities already known 
from the EOS determination. For example, the needed quantities are the chemical po- 
tentials, effective masses and single particle potentials of the baryon components. The 
first of these are generally tabulated in the EOS itself. In the nonrelativistic formalism, 
the effective masses and single particle potentials are simple functions of the particle 
number densities and kinetic energy densities. In the relativistic formalism, the effec- 
tive masses and single particle potentials are straightforwardly related to the particle 
number densities, internal energy and pressure by the minimization conditions. 

In the near future, we will prepare tables of both nonrelativistic potential and relativistic 
mean-field models of the EOS combined with their opacities, and make these generally 
available. Models including hyperons, kaon condensates and/or quarks will be computed. 

It is difficult to state categorically how the opacities reported here will affect the current 
generation of supernova or protoneutron star calculations. For example, although it seems 
clear that the addition of hyperons in matter will decrease mean free paths for a given 
density and temperature, the hyperon-bearing matter has a smaller temperature for a given 
entropy and density. Smaller temperatures tend to increase the mean free paths, and this 
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effect appears to more than cancel the original decrease. The feedback between opacities 
and EOS in a given astrophysical setting must be calculated consistently. 

It must be stressed that while we have included some effects of strong interactions through 
density dependent single particle excitations, there remain important collective effects from 
both density and spin/isospin dependent excitations |I0| , [I7P7| and from other density de- 



pendent in-medium correlations |pq| . The magnitude of these effects has so far only been 
assesed in some special cases such as nondegenerate symmetric matter or pure degenerate 
neutron matter or neutrino-poor beta-equilibrated matter. Furthermore, coupling to the 
A(1230) isobar |57| and screening by virtual particle-hole pairs created in the final state 
interactions |IB] could reduce the effective matrix elements. The net impact of such effects 
could be a reduction of some cross sections by as much as a factor of 2 at high densities. 
Work is in progress to calculate the role these effects on the neutrino cross sections at all 
temperatures of relevance and for all possible compositions, and will be reported subse- 
quently. 
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APPENDIX A: POTENTIAL MODELS 



Here, we outline a potential model for a system of neutrons and protons at finite tem- 
perature. With suitable choices of finite range interactions and parameters, this model re- 
produces the results of the more microscopic calculations (for more details, see, for example, 
Ref. |[47|). We begin with the energy density 
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where the factor 2 denotes the spin degeneracy and fi for i — n,p are the usual Fermi- 
Dirac distribution functions and m is the nucleon mass. It is common to employ local contact 
interactions to model the nuclear potential. Such forces lead to power law density-dependent 
terms in V(n). Including the effect of finite-range forces between nucleons, we parameterize 
the potential contribution as 
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where x = n p /n and u = n/n , with n denoting equilibrium nuclear matter density. 
The function g(k, A«) is suitably chosen to simulate finite range effects. The constants 
A, B,a, Ci, and C*2, which enter in the description of symmetric nuclear matter, and the 
additional constants xq, x$, Zi, and Z 2 , which determine the properties of asymmetric 
nuclear matter, are treated as parameters that are constrained by empirical knowledge. 

The single particle spectrum entering the Fermi-Dirac distribution functions fi may 
be written as 



21.2 



Ci(k) 



h z k 

2m,- 



+ Ui(k; n,x,T) , 



(104) 



where the single particle potential Ui(n,x,k;T), which is explicitly momentum dependent, 
is obtained by a functional differentiation of the potential energy density in Eq. (|103|) , with 
respect to the distribution functions /j. Explicitly, 
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where the upper (lower) sign in =p is for neutrons (protons) and i ^ j- 
The Landau effective mass is defined through the relation 
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where the prime denotes a derivative with respect to momentum. The finite range interac- 
tions may be approximated by effective local interactions by retaining only the quadratic 
momentum dependence: i.e., g(k,Ai) = 1 — (k/Ai) 2 . The energy density in Eq. (|103|) and the 
single particle potential in Eq. ( |105| ) then take the forms stemming from Skyrme's effective 
interactions ||34)| . Combining the finite range quadratic momentum term with the free ki- 
netic energy term, the single particle spectrum in Eq. ( J104 ) may be written for Skyrme-like 
interactions as 



ei(k) 



h 2 k 2 



2m 



- + Ui(n,x,;T) . 



(107) 



This resembles the free particle spectrum, but with a density dependent effective mass 
m* given by Eq. ([106]) with g'(k,Ai)\ kF . = l/( J R l 2 ^ 0) ), where Ri = Ki/{hkf) and 

EP = (hkf) 2 / (2m) is the Fermi energy of symmetric nuclear matter at the equilibrium 
density. Other forms of g{k, Aj), with more than a quadratic momentum dependence, also 
offer a viable description of the energy density and the single particle potential (see for 
example, Ref. ( [0)), but do not lead to a spectrum resembling the free particle spectrum 
in Eq. ( |107| ). This simple form of the spectrum allows a direct evaluation of the neutrino 
opacities including the effects of interactions along the lines developed for noninteracting 
baryons in Sec. III. For more general momentum dependent interactions, the evaluation of 
the neutrino opacities requires more complicated techniques which take into account the full 
momentum dependence in the four-momentum conserving delta functions. 

The parameters A, B,a, C±, and C2, are determined from constraints provided by the 
empirical properties of symmetric nuclear matter at the equilibrium density no = 0.16 fm~ 3 . 
With appropriate choices of the parameters, it is possible to parametrically vary the nuclear 
incompressibility K so that the dependence on the stiffness of the EOS may be explored. 
In this work, we have chosen K Q = 180 MeV for which 



A = 159.47 , B = -109.04 , a = 0.844 , G x = -41.28 and C 2 = 23 . 



(108) 
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Except for the dimensionless a, all quantities above are in MeV. The finite-range parameters 
A 1 = l.hpf and A 2 = 3pP . 

In the same vein, by suitably choosing the parameters x , x 3 , Zi, and Z 2 , it is possible 
to obtain different forms for the density dependence of the symmetry energy S(n) defined 
by the relation 



E(n,x) = e(n,x)/n = E(n, 1/2) + S{n)(l - 2x) 2 + 



(109) 



where E is the energy per particle, and x = n p /n is the proton fraction. Inasmuch as the 
density dependent terms associated with powers higher than (1 — 2a;) 2 are generally small, 
even down to x — 0, S(ri) adequately describes the properties of asymmetric matter. The 
need to explore different forms of S(n) stems from the uncertain behavior at high density 
and has been amply detailed in earlier publications |]45 , |48 | 
the potential part of the symmetry energy to vary as u. For this case, 



In this work, we have chosen 



x = -0.410, x 3 = -0.5 Z x = -11.56 MeV, and Z 2 = -4.421 MeV 



(110) 



Since repulsive contributions that vary faster than linearly give rise to acausal behavior 

This may be 



at high densities, care must be taken to screen such repulsive interactions [4S 



achieved by dividing the term proportional to u a+1 (when o > 1) by the factor 



! + §£' 



I" (* + x 3 ){l-2x) 



u 



er-l 



111) 



where B' is a small parameter introduced to maintain causality. Note that the single 
particle potential in Eq. ( |105| ) must then be accordingly modified. The appropriate terms 
may be obtained by using the relations (d/dn n )\ np = d/dn — (x/n)d/dx for neutrons and 
(d/dn p )\ nn = d/dn + [(1 — x)/n]d/dx for protons. 

For a fixed baryon density n, proton fraction x, and temperature T, an iterative procedure 
may now be employed to calculate the density dependent single particle potentials U n and U p , 
and the chemical potentials fi n and \i v . The calculational procedure is detailed in Ref. p? 
These quantities, in conjunction with the requirements of chemical equilibrium in Eq. (| 



and charge neutrality in Eq. ( |j"U| ) determines the composition of stellar matter at finite 
temperature. The specification of the spectrum in Eq. ( |107| ) then allows for a calculation of 
the neutrino opacities which are consistent with the underlying EOS. 
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APPENDIX B: EFFECTIVE FIELD-THEORETICAL MODELS 



In a Walecka-type relativistic field-theoretical model the interactions between baryons 
are mediated by the exchange of a, lu, and p mesons. The Lagrangian density is given by |35| , 



L — Lh + Lt 

= ^Bi-iYdf, - g^Bl^^p, - gpBrfhp • t - M B + g<,B(y)B 

B 

+ l -d,ad»o - \m\a 2 - U{a) 

+ Y,l{-i^d li -m l )l. 

i 

Here, B are the Dirac spinors for baryons and t is the isospin operator. The sums include 
baryons B = n,p, A, £, and S, and leptons, I = e~ and p~ . The field strength tensors for 
the co and p mesons are W pv = d^uj v — d v uj^ and B^ = <9 M b„ — d u b p , respectively. The 
potential U(a) represents the self-interactions of the scalar field and is taken to be of the 
form 

U(a) = ^bM n ( 9(7N a) 3 + ^c(g aN a) 4 . (112) 

Electrons and muons are included in the model as noninteracting particles, since their in- 
teractions give small contributions compared to those of their free Fermi gas parts. 

In the mean field approximation, the partition function (denoted by Zh) for the hadronic 
degrees of freedom is given by 

lnZ H = (3v[\mlul + \m 2 p bl-\mla 2 -U(o) 



B J 



|£ In (l +»-**.->), (113) 



where (3 = (kT) 1 and V is the volume. The contribution of antibaryons is not significant 
for the thermodynamics of interest here, and is therefore not included in Eq. ( |113| ). Here, the 



effective baryon masses M B = Mb — g a B& and E* B = yk 2 + M B 2 . The chemical potentials 
are given by 

Pb = v b + guBU + gpBkabo , (114) 

where t 3B is the third component of isospin for the baryon. Note that particles with t 3 g = 0, 
such as the A and E° do not couple to the p. The effective chemical potential v B sets the 
scale of the temperature dependence of the thermodynamical functions. 

Using Zh, the thermodynamic quantities can be obtained in the standard way. The 
pressure P# = TV^ln Z H , the number density for species B, and the energy density e# 
are given by 
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n B 



d3k ( e PiE* B -u B ) , l Y 
(2tt)3 V e + L ) 



en = \mW + U(a) + \m^l + fm^ + 2^ / (S E ^ + if' . (115) 

The entropy density is then given by sh = P(£h + Ph — J2b Hb^b)- 

The meson fields are obtained by extremization of the partition function, which yields 
the equations 



m 2 . 



u o = Yl 9ojBn B ; m 2 p b = ^ g pB hBn B 

B B 

dU ^) , v- f d 3 k M* B fM _ u s , lV -i 



"V = — + Eto (^W + !)-' . (U6) 

The total partition function Z to4a ^ = ZhZ^, where is the standard noninteracting parti- 
tion function of the leptons. 

The additional conditions needed to obtain a solution are provided by the charge neu- 
trality requirement, and, when neutrinos are not trapped, the set of equilibrium chemical 
potential relations required by the general condition 

Hi = bifJL n - qifMi , (117) 

where b{ is the baryon number of particle i and qi is its charge. For example, when i = e~ , 
this implies the equalities 

Ha = Hyp = He = Hn > 
Hy- = Hz- = Hn + He, 

H P = Hy+ = Hn - He- (118) 



In the case that the neutrinos are trapped, Eq. (|117|) is replaced by 



Hi = biHn - qi(Hi ~ Hu e ) ■ (119) 

The new equalities are then obtained by the replacement \x e — > \x e — fi Ue in Eq. ( |118| ). The 
introduction of additional variables, the neutrino chemical potentials, requires additional 
constraints, which we supply by fixing the lepton fractions, Yli-, appropriate for conditions 
prevailing in the evolution of the protoneutron star. The contribution to pressure from 
neutrinos of a given species is P u = (l/247r 2 )/^. 

In the nucleon sector, the constants g a N, g^N, g P N, b, and c are determined by reproducing 
the nuclear matter equilibrium density tiq = 0.16 fm -3 , and the binding energy per nucleon 
(~ 16 MeV), the symmetry energy (~ 30 — 35 MeV), the compression modulus (200 MeV < 
K < 300 MeV), and the nucleon Dirac effective mass M* = (0.6 - 0.7) x 939 MeV at n . 
Numerical values of the coupling constants so chosen are: 

gaN/m a = 3.434 fm, g uN /m u = 2.674 fm, g pN /m p = 2.1 fm, 

6 = 0.00295, and c = -0.00107. (120) 
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These couplings yield a symmetry energy of 32.5 MeV, a compression modulus of 300 MeV, 
and M* /M = 0.7. This particular choice of model parameters are from Glendenning and 
Moszkowski |58[ and will be referred to as GM1 hereafter. The prevalent uncertainty in the 
nuclear matter compression modulus and the effective mass M* does not allow for a unique 
choice of these coupling constants. The high density behaviour of the EOS is sensitive to 
the strength of the meson coupling constants employed. Lacking definitive experimental and 
theoretical constraints, this choice of parameters may be considered typical. 

The hyperon coupling constants may be determined by reproducing the binding energy 
of the A hyperon in nuclear matter |]58|| . Parameterizing the hyperon-meson couplings in 
terms of nucleon-meson couplings through 

x aH = 9uHl 'daN, X w n = 9wH I ' 9uN ', % pH = 9pIl/9pN > (121) 

the A binding energy at nuclear density is given by 

(B/A) A = -28 = x uA g uN uo - x aA g aN cr , (122) 

in units of MeV. Thus, a particular choice of x a \ determines x u a uniquely. To keep the 
number of parameters small, the coupling constant ratios for all the different hyperons are 
assumed to be the same. That is 

X a = X aA = X a Y, = X a ~ = 0.6 , (123) 

and similarly for the u 

Xu, = x w a = Zu>e = x^e = 0.653 . (124) 
The p-coupling is of less consequence and is taken to be of similar order, 
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TABLE I: Charged current vector and axial couplings. Numerical values are quoted using 
D=0.756 , F=0.477, sin 2 9 W =0.23 and sin 2 # c = 0.053 (see Ref. g|). The couplings for the 
same reactions involving antineutrinos are identical, and i stands for e— , \x— or r— type 
leptons. For corrections arising due to explicit SU(3) breaking terms, see Ref. p0| . 



Reaction 




9v 


9a 


AS 


V(_ + n — ► 




1 


D + r = 1.23 





vt + £~ - 


>r + a 





^2/3D = 0.62 





V£ + E~ - 


> r + s° 


V2 


y/2F = 0.67 





v t + S° -» 


r + e+ 


-V2 


-v^F = -0.67 





z/^ + A -> 


r + e+ 





-^2/3D = -0.62 





z/^ + A -> 




-^3/2 


/3/2(F + D/3) = 0. 


39 1 


^ + E° -> 


r + p 


1 


^1/2D = 0.54 


1 


^ + £~ - 


> +n 


-1 


D-F = 0.28 


1 


^ e + /i _ — 




1 


1 






TABLE II: Neutral current vector and axial couplings. Neutral current couplings with 
baryons of all neutrino species, including antineutrinos, are identical, and £ stands for 
e— , or r— type neutrinos. Neutrino interactions with leptons have the same matrix 
elements as those with antineutrinos of the same flavor. For corrections arising due to 
explicit SU (3) breaking terms, see Ref. [[4(J . 

Reaction cy ca 

v~ e + e~ -> v~ e + e~ 1 + Asm 2 6 w = 1.92 1 

v^ + \l~ -> + \l~ 1 + 4 sin 2 6 W = 1.92 1 

v e + \l~ -> v e + fT -1 + Asm 2 8 w = -0.08 -1 

"h,t + e" -> ^ M ,r + e~ -1 + 4 sin 2 6» w = -0.08 -1 

v £ + n^v e + n -1 -D — F = -1.23 

z^+p^^+p 1 -4sin 2 ^ = 0.08 F + F = 1.23 

z/ £ + A^z/ £ + A -1 — F — D/3 — —0.73 

z/ £ + S-^z/ £ + S- -3 + 4 sin 2 % = -2.08 D — 3F = —0.68 

z/ £ + £+ -> ^ + £+ 1 -4sin 2 # w = 0.08 D+F =1.23 

v e + E° -> v e + E° -1 D- F = 0.28 

14 +- S~ z/< +- S~ -3 + 4 sin 2 ^ = -2.08 D - 3F = -0.68 

z/ £ + S°^^ + S° -1 -D-F = -1.23 

^ + E° -»• v t + A 2^/^ = 0.87 
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FIGURE CAPTIONS 



Fig. 1. Lowest order Feynman diagram for v% + — > I + B4. The symbols B{ and I denote 
baryons and leptons, respectively. P, are the particles' four- momenta and g M = (qo,q) is the 
four- momentum transfer. Fig. (la) is the absorption reaction and Fig. (lb) is the scattering 
reaction. 

Fig. 2. Composition and chemical potentials for noninteracting matter in beta equilibrium 
at Yl = 0.4 at different temperatures. Top panels: Individual concentrations Y> = rii/ns, 
where i = n, p, e~ and v e . Bottom panels: The lepton chemical potentials and \i = ji n —ji p = 

Fig. 3. Absorption mean free paths A = a/V in noninteracting matter in beta equilibrium 
at different temperatures. The neutrino energy is taken to be the Fermi energy (E v = fj, u ). 
The solid curves show exact results from Eqs. (|T5|) and (|2T|), the dashed curves are from the 
degenerate approximation (Eq. (|26f) ), the dot-dashed curves are the nondegenerate results 
from Eq. (|3"TD, and the long dashed curves represent the elastic approximation from Eq. (pip. 

Fig. 4. Comparison between scattering and absorption reactions at T = 10 MeV for 
neutrino-trapped matter (Yl = 0.4), assuming E u = n v . The neutrino energy was set equal 
to the local neutrino chemical potential. In the left panels the contribution to the neutrino 
mean free path due to individual reactions are shown. In the right panel, the solid curve is 
the ratio of the scattering to absorption mean free paths and the the dashed curve is the 
neutron scattering/absorption ratio. 

Fig. 5. Compositions and chemical potentials for noninteracting matter in beta equilibrium 
with no trapped neutrinos (Y v = 0) at different temperatures. Top panels: Individual 
concentrations Y^ = rii/nB, where % = n, p, and e~. Bottom panels: The electron chemical 
potential \i e = fi n — fi p . 

Fig. 6. Absorption mean free paths A = a/V in noninteracting matter in beta equilibrium 
with no trapped neutrinos at different temperatures. The neutrino energy is taken to be 
Ei = 3T. The leftmost panel clearly shows that at low temperatures the mean free path is 
very large since the reaction is kinematically supressed. 

Fig. 7. Temperature dependance of the neutrino mean free path at nuclear saturation 
density for Y v = and E v = 3T. Left panel: Individual contributions due to scattering and 
absorption reactions to the mean free path, the long-dashed curve is the total scattering 
opacity. Right panel: Ratio of the total scattering opacity to that of absorption. 

Fig. 8. Particle fractions (upper panels) and lepton chemical potentials (lower panels) in 
the nonrelativistic potential Skyrme model for Yl = 0.4 at various temperatures. Note that 
the proton fraction is larger than in the noninteracting case. 

Fig. 9. Neutrino absorption mean free path in neutrino-trapped matter with Yl = 0.4. 
Results for the limiting cases discussed in the text are also shown, with the same notation 
as in Fig. 3. 
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Fig. 10. Neutron (solid curve) and proton (dashed curve) effective masses as a function of 
baryon density in the Skyrme model. The curves shown are for T = 30MeV. 

Fig. 11. Compositions (upper panels) and electron chemical potentials (lower panels) for 
neutrino-free matter at various temperatures in the Skryme model. Note that the proton 
fraction is larger when compared with the noninteracting case and, at low densities, is quite 
sensitive to the temperature. 

Fig. 12. Upper panel: Neutrino absorption mean free path in neutrino-free matter with the 
Skyrme model, assuming E v = 3T. The exact and limiting cases discussed in the text are 
shown, with the same notation as in Fig. 3. For T = 5 MeV, the elastic approximation is 
off-scale. Lower panel: Comparison of absorption and scattering mean free paths for this 
case. 

Fig. 13. Compositions (upper panels) and chemical potentials (lower panels) in field- 
theoretical model for nucleonic matter in beta equilibrium with fixed lepton number 
Yl = 0.4. 

Fig. 14. Absorption (upper panels) and scattering (lower panels) mean free paths A = a/V 
for neutrino energy E u = {i Ue in beta equilibrium matter with Yl = 0.4, for the field- 
theoretical model. 

Fig. 15. Composition (upper panels) and chemical potentials (lower panels) matter in beta 
equilibrium with no trapped neutrinos (Y u = 0) for the field-theoretical model. 

Fig. 16. Absorption (upper panels) and scattering (lower panels) mean free paths A = a/V 
for neutrino energy E u = 3T in matter with no trapped neutrinos for the field-theoretical 
model. 

Fig. 17. Comparision with previous cross sections employed in protoneutron star calcula- 
tions. The upper panels show comparisions with limiting forms valid for the nondegenerate 
(left) and degenerate (right) matter. The lower panels show comparisions with interpolated 
cross sections; the interpolation scheme in Eq. fl8"6"|)is shown in the left panel and that in 
Eq. flSTD is shown in the right panel. 

Fig. 18. Compositions (upper panels) and chemical potentials (lower panels) for neutrino- 
trapped matter (Yl = 0.4) including hyperons in the field-theoretical model. // = \i n — \i p = 

Fig. 19. Absorption mean free paths A = (a/V)^ 1 for neutrino energy E v = \i Ve in matter 
with fixed lepton number (Yl = 0.4) in the field-theoretical model. Upper panels show both 
the results for matter without hyperons (dashed curves) and with hyperons (solid curves). 
Lower panels show the relative contributions to the absorption mean free path for matter 
with hyperons. 

Fig. 20. Scattering mean free paths for neutrino energy E v = fi Ue in matter with fixed lepton 
number (Yl = 0.4) in the field-theoretical model. Upper panels show the results for the toal 
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scattering mean free path for matter without hyperons (dashed curves) and with hyperons 
(solid curves). Lower panels show the relative contributions to the scattering mean free path 
for matter with hyperons. 

Fig. 21. Compositions (upper panels) and chemical potentials (lower panels) for neutrino- 
free matter (Y u = 0) including hyperons in the field-theoretical model. 

Fig. 22. Absorption mean free paths A = (a/V)^ 1 for neutrino energy E u = 3T in matter 
with no trapped neutrinos (Y v = 0) in the field-theoretical model. Upper panels show both 
the results for matter without hyperons (dashed curves) and with hyperons (sold curves). 
Lower panels show the relative contributions to the absorption mean free path for matter 
with hyperons. 

Fig. 23. Scattering mean free path for neutrino energy E v = 3T in matter with no trapped 
neutrinos (Y u = 0) in the field-theoretical model. Upper panels show the results for the total 
scattering mean free path for matter without hyperons (dashed curves) and with hyperons 
(solid curves). Lower panels show the relative contributions to the scattering mean free path 
for matter with hyperons. 

Fig. 24. Neutrino diffusion coefficient D 2 defined in Eq. ( |9lD in matter with and with out 
hyperons for Yl = 0.4 in the field-theoretical model. The left panel shows D 2 as a function 
of baryon density in matter containing only nucleons and leptons. The ratio of D 2 in matter 
without hyperons to D 2 in matter with hyperons is shown in the right panel. 

Fig. 25. Neutrino diffusion coefficient -D4 defined in Eq. fl9T|) in matter with and with out 
hyperons for Y v = in the field-theoretical model. The left panel shows -D4 as a function of 
baryon density in matter containing only nucleons and leptons. The ratio of D4 in matter 
without hyperons to -D4 in matter with hyperons is shown in the right panel. 
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